Trimer liquids and crystals of polar molecules in coupled wires 
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We investigate pairing and crystalline instabilities of bosonic and fermionic polar molecules con- 
fined to a ladder geometry. Combining analytical and numerical techniques, we show that gases of 
composite molecular dimers as well as trimers can be stabilized as a function of the density difference 
between the wires. A shallow optical lattice can pin both liquids, realizing crystals of composite 
bosons and fermions. We show that these exotic quantum phases are robust against conditions of 
confinement of the molecular gas to harmonic finite-size potentials. 
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Binding among three or more particles plays a fun- 
damental role in various physical systems ranging from 
quark confinement in QCD to tetramer bi-exciton for- 
mation in carbon nanotubes [l[; in quantum magnets it 
is responsible for the formation of many-body quantum 
phases of spin-multipoles 0. In atomic quantum gases, 
two-body pairing is at the heart of paradigmatic phenom- 
ena, such as the observation of the BEC-BCS crossover Q 
in mixtures of ultracold fermionic atoms and the predic- 
tion of FFLO Q and Sarma Q phases in spin- unbalanced 
gases. The study of few-body pairing^ however, is usually 
confined to Efimov-like resonances [6| , while the realiza- 
tion of finite-density liquids made of composite particles 
is hindered by losses due to three-body recombination. 
Polar molecules confined to low-dimensional geometry 
provide a new opportunity to study inter-molecular pair- 
ing mechanisms and the associated quantum phases in 
a setup where collisional losses and also chemical reac- 
tions arc suppressed Pairing of two spin-polarized 
fermionic molecules across coupled two-dimensional (2D) 
layers [13] or one-dimensional (ID) wires 11 1 has already 
lead to the prediction of, e.g., 2D inter-layer superfluid- 
ity for the special case where the number of molecules 
is the same in all layers (wires). However, the pairing 
dynamics in the general situation where this number can 
vary across the layers (wires), as it happens in experi- 
ments has so far remained largely unexplored. Dif- 
ferent particle numbers in the layers lower the overall 
(spin-rotational) symmetry of the problem, which makes 
the underlying physics drastically different [13l | : a general 
pairing mechanism may exist for the formation of stable 
multi-molecule composite structures. The goal is now 
to determine whether (spin-rotational) symmetry break- 
ing induces stable multimer liquids in confined quantum 
gases, where few-body processes are usually associated 
with resonances and losses 

In this letter we study a situation where polar 
molecules are confined to two coupled ID wires, under 
conditions where the number of molecules can be the 
same or vary across the wires, similar to current experi- 
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FIG. 1: (color online), (a): Molecules with dipole strength 
d a (a — 1,2) in a two-leg ladder with interwire distance g 
and dipolar interactions V(r) = di<fe(l — 3cos 2 #)/r 3 . (b- 
d): Sketch of phases (see text): (b) dimer crystal in a lattice 
with spacing A/2 and hopping energy t a ; (c) dimer liquid; 
(d) trimer liquid, (e-f): Numerical results for density distri- 
butions n a vs position x in the presence of an optical lat- 
tice along the wires [see panel (b)] and weak harmonic con- 
finement, showing trimers. Parameters [see also Eq. ((2} and 
text]: g = A/2, hopping asymmetry £2/^1 = 0.01, di/do = 3, 
d-ifda = 0.5 with d\ — tig 3 ; the curvature of the harmonic 
potential is /C/ti = 0.03 (see text). Particle numbers N a : (e) 
N 2 = 2Ni = 12; (f) N 2 = 2JVi + 2 = 14; as a reference, blue 
points denote 2ni. The trap center is at x/(X/2) — 60. 



ments. This presents several new scenarios: (i) when the 
population is identical in the wires, a two-body bound- 
state is always present and is responsible for the appear- 
ance of dimer liquids, similar to 2D [l2|, see Fig. QJc). 
However, we find that (ii) few-body pairing is favored for 
any ratio of populations between the wires P = p/q, with 
p,?eN, which is a dense set between ]0, 1]. This leads to 
the stabilization of novel quantum phases of interacting 
fermionic or bosonic composite particles, where repulsive 
intra- wire interactions ensure collisional stability. In par- 
ticular, we prove that a gas of trimers made of two par- 
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tides on one wire and one in the other may be stabilized 
in these systems, Fig. [ljd). Both dimer and trimer liq- 
uids may be pinned by a weak period potential commen- 
surate with the density, leading to a Luttinger-staircasc 
of quasi-lD composite crystals Q. By numerical simu- 
lations, we show that these exotic quantum phases are 
robust to conditions of trapping in parabolic finite-size 
potentials and number fluctuations in the wires: unusual 
" wedding-cake" structures are formed where heavy com- 
posite particles are flanked by lighter one, Fig. QJe-f). 

Our starting point is the Hamiltonian H = Ha + 
H12 for the dynamics of molecules in the configuration of 
Fig. [TJ with a = 1, 2 and H a the single- wire term 



H a = J dx tpi(x) 
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and H\2 = (didz) / (8n) J dx dx'V{x — x')n\{x)ri2{x') the 
inter- wire coupling, with V(x — x') = [1 — 3 cos 2 (6)]/[g 2 + 
(x — x') 2 ] 3 / 2 showing a short-distance inter-wire attrac- 
tion. Here, m a ,d a are the mass and the dipole strength 
respectively, and ip a (4>a) are fermionic or bosonic anni- 
hilation (creation) operators, with n a (x) = ipl l (x)4'a(x); 
U a (x) = U a sm 2 (2Trx/\) represents an underlying pe- 
riodic potential, as usually provided by an optical lat- 



tice [Toj with wavelength A and depth U a \ 6 is the an- 
gle between particles in different wires, with distance g, 
Fig. [TJ The ID Hamiltonian H is valid for interparticle 
distances n" 1 3> (d 2 1 /huj±) 1 ^ 3 with uj± the frequency of 
transverse confinement provided by, e.g., a 2D optical lat- 
tice. For LiCs, RbCs and KRb molecules with d\ = 5.6, 
1.25, and 0.5 Debye, respectively, and uj± = 2irx 100 kHz, 
(dl/huj^) 1 ^ is of the order of 360, 130 and 70 nm 0. 
In addition, <ii<i 2 /<? 3 < U a , [lf|. Similar setups may be 
obtained on chip-based microtraps (lij |. 

In it is shown that in the absence of an optical lat- 
tice [U a = 0) and of inter-wire interactions, the dynam- 
ics in each wire is described by an effective Tomonaga- 
Luttinger liquid (TLL) theory with Hamiltonian [l[ 

n ° = ~2~n J dX [( d M x V 2 / K <* + K*(dMx)) 2 ] ■ 

Here v a and K a = (1 + O^n^Ac,) -1 / 2 arc the ef- 
fective sound velocity and the TLL parameter, respec- 
tively with R a — m a d 2 a /(2ttH 2 ) the intra-wirc dipole 
length, and i? Q , <p a represent long- wavelength density and 
phase fluctuations [1|. For finite interactions, Eq. is 
then effectively described by % — J2a^ a + ^ 12 j w ith 
H12 = Hf+Hb- Here, in the weak coupling regime we ob- 
tain Hf = — 24gVii / dxd x '&i(x)d x '&2{x) for the quadratic 
forward-scattering part of the interactions [l| , by approx- 
imating the interwire interaction with its zero-component 
Fourier transform[l|; Hb is the back scattering part, with 
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FIG. 2: (color online), (a): Phase diagram in the balanced 
case with d a = d, t% = ii and n a = 0.2 on a lattice with 
spacing A/2 and d%> = ii(A/2) 3 . Green triangles, red circles 
and blue squares are numerical results for dimer liquid (DL), 
charge-density wave (CDW) and dimer crystal (DC) phases, 
respectively; dashed and dashed-dotted lines are qualitative 
phase boundaries. Inset: scaling of the dimer-pairing spin- 
gap A s as a function of d 2 for several example densities n a 
and interwire distances g: [n a = 0.2, g = 0.35A] (red points), 
[0.4, A/2] (green triangles) and [0.2, A/2] (black diamonds), 
(b): Dimer correlation function T>[x) vs x in the phases of 
panel (a), with g/(X/2) = 0.5, n a — 0.2. Top to bottom (con- 
tinuous lines): d/d = 0.5 (DL), 1.3 (CDW) and 2.2 (DC); the 
dashed-dotted line is the CDW-DC quantum phase transition 
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a typical sine-Gordon-type (sG) form, to be discussed be- 
low. The effect oiHf is to modify the effective TLL pa- 
rameters, while Wb can induce novel pairing instabilities. 
Here, we first discuss the case of two identical [balanced) 
coupled wires, and then the more general case of [unbal- 
anced) wires with different densities, molecular masses 
and interactions. 

In the balanced case, the quadratic part ^2 a T-La +Hf 
of T-L can be diagonalized by introducin g st andard charge 
and spin fields, i? CiS = ($i ± $2) /a/2 [20(, leading to a 
description in terms of coupled TLLs with effective pa- 
rameters K c s . For weak interactions we estimate K CiS = 

K x [1 TFi 2 A'i/ui]~ 1/2 , L12 = dl/[2An 2 hg 2 ). The term 

Hb has the sG form Hb oc — / dx cos[2^/2d s [x)], 

and, in agreement with Berezinskii-Kosterlitz-Thouless 
(BKT) theory [J|2l[, it is relevant, thus opening a spin 
gap (e.g., pairing), if K s < 1; this condition is always 
satisfied and shows that, similar to the 2D case, pairing 
of molecules across the wires is always favored in ID, 
even for an infinitcsimally small attraction between the 
wires with g, R a <C n^ 1 . The role of dipolar interactions 
is evident in the charge sector: in contrast to models 
with attractive contact interactions such as the Hubbard 
model [l[, here K c can be much smaller than 1; as an 
example, in the strongly interacting regime n a R a 3> 1, 
where dimers are well approximated by tightly-bound 
composite particles with effective mass M = 2m and 
dipole strength D = 2d, K c ~ 2/(1 + 5.8n Q i? Q ) 1/2 . As 
a result, the many-body groundstate shows a crossover 
from a dimer liquid (DL) with dominant pair corre- 
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lations V(x) = {tpl^l^ipi^+x^i+x) ~ \x\ 1/Kc for 
K c > 1 to a charge-density wave (CDW) with domi- 
nant density correlations Q(x) = (mni+ x ) ~ l^l - ^ for 
K c < l,m = rii t i + n,,2- In this regime, the CDW 
can be pinned by a very shallow optical lattice commen- 
surate with the particle density, stabilizing a Luttinger 
staircase of dimer crystals (see [20[ for more details). 

We verify numerically these predictions in the deep 
lattice regime U/E r ^> 1, with E r the lattice recoil, 
where an appropriate description is given in terms of an 
anisotropic extended Hubbard model (AEHM) [j^ : 

x - ^ / t \ t 2d\d 2 v~^v 
H = - y J a {c l a i c a ,j + i + h.c.) p— y j ni.in 2 ,i + 
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which we analyze using a quasi exact Density-Matrix- 
Renormalization-Group technique (DMRG) [23j. Here, 
Vij describes the anisotropic part of the dipolar interac- 
tion, t a = 1 sets the energy scale. In the balanced case, 
d a = d and n a — n. The field theoretical description of 
the AEHM in terms of continuum fields is equivalent to 
the one in the limit of a shallow lattice given above, and 
we thus expect a similar qualitative behavior. 

Figure [5] shows the phase diagram for a commensurate 
density n a = 0.2. By fixing g = A/4 and increasing d, 
we find first a crossover from a TLL of dimcrs (DL) to 
a CDW, and then a BKT-type pinning quantum phase 
transition to a DC with n a = 1/5 (phase boundaries are 
discussed in [2(|). Examples of V(x) are plotted in panel 
(b) for all three phases, where the dash-dotted line marks 
the transition from power-law to exponential decay. We 
also calculate the spin gap A s , by performing a finite- 
size scaling of A S (L) = E L (N,N) - E L (N + 1,N- 1), 
with El(M,M') the ground state energy at finite size 
L in the sector with n\ — M,ri2 = M', for different 
densities and g [Inset of panel (a)] . We find that a finite 
gap is present in the entire phase diagram, as expected, 
although it is small for weak interactions due to the BKT 
scaling A s cx exp[— (3d 2 }, [l[. 

The unbalanced case presents unconventional instabil- 
ities. As shown in [l3j], the Haldane expansion of density 
operators in a two-component TLL generates a infinite 
series of massive terms coming from different combina- 
tions of vertex operators 



%b = J]] Gp q J dx cos[2x(pfci?i 



-qk F2 ) + 2(p9 1 -q9 2 )}, 



where G PtQ are model dependent coefficients and kp a = 
7rn a . With the exception of the simplest case p = q = 1, 
these terms are usually negligible from a renormalization 
group point of view. However, we find that strong dipo- 
lar interactions drastically enhance the effect of % [20|, 
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FIG. 3: (color online), (a): Phase diagram in the unbal- 
anced case, as a function of d\/a% and £2/^1. Here, di = 6^2, 
g = A/2, n 2 = 2m = 0.1 and do = £i(A/2) 3 ; circles, tri- 
angles and squares are numerical results for the trimer liquid 
(TL), phase separation (PS) and two independent Tomonaga- 
Luttinger-Liquid phases (2TLL), respectively; lines are guides 
for the eye. Inset: Trimer binding energy As/ii vs £2/^1 for 
several choices of (di/do; ife/do): triangles, squares and circles 
are (0.5; 3), (0.3; 1.8) and (0.5, 2), respectively, (b-c): Corre- 
lation functions for trimers T{x) (dashed-dotted line), inde- 
pendent TLL Ci(x) (continuous line) and dimer liquid T>(x) 
(dashed line) vs 1 in a system of size L — 120 (see text), (b): 
Correlations in a TL with di/do = 3, d2/do = 0.5, £2/^1 = 0.1; 
(c) 2TLL with di/do = 2.4, d 2 /d = 0.4, t 2 /ti = 0.4. 



allowing the formation of multiparticle composites [13| . 
or multimers, made of p/q particles on the upper/lower 
wire respectively. In particular, composite objects with 
p = 1 may become particularly stable, corresponding to 
a population ratio n 2 /n\ = k € N; in such a situation, 
analogous to the above discussion of dimers, a term in Hi, 
with the sG form G Kj i J dx cos[2(ki?i — $2)] mav become 
relevant and stabilize a multimer liquid, uniquely identi- 
fied by the finite gap associated with the bound state for- 
mation and an algebraic decay of multimer correlations 
(M^(0)M(x)), M = (V>i) K -i/> 2 , while the single-particle 
and dimer correlations T>(x) decay exponentially. How- 
ever, qualitative estimates (see [20|) indicate that com- 
paratively large interaction strengths are needed in order 
to stabilize such a liquid, and the critical strength in- 
creases with increasing k. Thus, below we focus on k = 1 
and we investigate numerically the possibility to realize a 
trimer liquid (TL) by considering the AEHM with both 
interaction and hopping asymmetry, quantified by the 
ratios d 2 jd\ and t 2 /t\. 

To get a guidance on possible trimer instabilities, we 
first compute the binding energy of the three body prob- 
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km, A B = \im L ^ oo [E L (l,l)+E L (0A)-E L (l,2)}. Based 
on a Born-Oppcnhcimcr approach [24J, we expect a size- 
able binding in the regime t\ ^> t 2 and d\ >d,2, where a 
fast particle in wire 1 binds two repulsive heavy molecules 
in wire 2. Numerical examples of A# vs t 2 /t\ are given in 
the Inset of Fig. for a few values of d 2 , d\ and g = X/2, 
as trimers are unfavored for too large or too small in- 
terwire distances [2(|. In particular, for strong attrac- 
tion with j « A, the formation of a dimer plus an un- 
paired particle is favored (see also the discussion of phase- 
separation (PS) below). We have further investigated the 
phase diagram in the low-density limit (n 2 = 2n\ = 0.2 
and 0.1) as a function of di, d 2 for system sizes up to 
L = 120 using DMRG. The TL phase is characterized by 
an exponential decay of both V(x) and single particle cor- 
relation function C a (x) = (c^ ^a ^), whereas the trimer 

correlator T(x) = (cl^.i^.i+iCa.i+s+i^.i+xCi.i+x) de- 
cays algebraically. This is in contrast to the case of two 
coupled TLL, where no correlation is exponentially sup- 
pressed. The phase diagram for d\ = 6d 2 ,n 2 = 0.1 is 
plotted in Fig. [3] (see also [20[ for the dependence on 
d\/d 2 ): the TL extends in a broad region, and survives 
even for comparatively small interaction strength and in- 
teraction asymmetry, albeit in both cases small hopping 
rates t 2 < 0.2 arc needed. Sample correlation functions 
for both TL and 2TLL arc plotted in FigEtb-c). For 
large interactions did 2 /g 3 3> t a , microscopic phase sep- 
aration (PS) can occur. The latter corresponds to the 
formation of a gas of strongly-bound, mutually repulsive 
dimers across the wires, coexisting with a gas of unpaired 
molecules in wire 2, see |20| . 

Effect of the trap. We have investigated the fate of 
these exotic trimer bound states in the lattice under con- 
ditions of harmonic trapping in finite-size potentials with 
Hamiltonian H t = (JC/ L 2 )Y Jl n i {L/2 - if, with JC the 
curvature. We characterize the TL phase in inhomo- 
geneous trapped systems via the local order parameter 
8iii = 2ni i — n 2 ^, describing the deviations between the 
density of the majority and minority components poj ]. 
Trimers exist if drii <C 1. We find that trimers are re- 
markably robust to the presence of high-density, finite- 
size effects, and fluctuations in particle numbers in the 
wires. The general situation is one where composite par- 
ticles with a larger mass (e.g., trimers) occupy the central 
region of the trap, and are flanked by lighter particles, be- 
ing dimers or single excess particles. This determines an 
unusual "wedding-cake" structure in the density profile, 
reminiscent of that observed with cold atom with contact 
interactions [L5J. Figures H{e-f) show example results of 



the density profile in the trap, in the presence of compos- 
ite trimers for of an exactly-matching number of 
particles in the two wires, and a situation where the par- 
ticle number in wire 2 is larger by 20 percent. In both 
cases, we find a well defined region at the trap center 
where trimers are formed, with Srii ~ 10 -8 . 



Significant asymmetry in dipoles and mass is obtained 
in mixtures of, e.g., LiCs and RbCs molecules trapped 
in independent optical lattices, where a weak optical 
lattice in the wire direction provides additional tun- 
ing of the (effective) mass. Alternatively, single-species 
molecules can be prepared in different internal (e.g., ro- 
tational, vibrational) states with different dipole mo- 
ments [!, [25|. For a ID lattice in the wire direction, 
the effective mass is tuned via internal-state-dependent 
tensor shifts 25, [26|. All phases discussed above can 



be detected via the measurement of decay of correlation 
functions [see, i.e., Fig. [3fb-c)1, using, e.g., direct in situ 
imaging techniques |27J]. Particle correlations across the 
wires will be detected via in situ imaging (27| as well as 
noise correlation measurements [l5|. In addition, the TL, 
DL and DC may be spectroscopically probed (l3|. For ex- 
ample, molecular dimers with sizeable, spectroscopically 
resolvable, gaps A s > 1/iK are obtained by trapping, e.g., 
LiCs molecules in a deep 2D optical lattice with spacing 
A/2 = 400, trapping frequency uj± = 2tt x 100 kHz. 

The realization of the above scenario in polar molecule 
experiments will lead to the first observation of trimers 
in the many-body context with cold gases. 

We thank E. Burovski and A. Chotia for discussions, E. 
Ercolessi and F. Ortolani for discussions and for provid- 
ing the DMRG code. This work was supported by MURI, 
AFOSR, EOARD, IQOQI, the Austrian FWF, the EU 
through NAME-QUAM, COHERENCE and AQUTE. 

Note: after completion of this work, we became aware 
of the related work [28| on few-body boundstates of polar 
molecules in coupled tubes. 



[I] R. Matsunaga, K. Matsuda and Y. Kanemitsu, Phys. Rev. 
Lett. 106, 037404 (2011). 

[2] J. Sudan et al, Phys. Rev. B 80, 140402 (2009) (R). 

[3] M. Inguscio, W. Ketterle, and C. Salomon (Eds.), Ultra- 
cold Fermi Gases, IOS Press, Amsterdam (2008). 

[4] P. Fulde and R. A. Ferrell, Phys. Rev. 135, A550 (1964); 
A. I. Larkin and Yu. N. Ovchinnikov, Zh. Eksp. Teor. Fiz. 
47, 1136 (1964). 

[5] G. Sarma, J. Phys. Chem. Solids 24, 1029 (1963). 

[6] F. Ferlaino and R. Grimm, Physics 3, 9 (2010). 

[7] H.P. Biichler et al, Phys. Rev. Lett. 98, 60404 (2007); A. 
Micheli et al, . ibid. 105, 073202 (2010); G. Quemener 
and J. L. Bohn, Phys. Rev. A 81, 022702 (R). 

[8] M. H. G. de Miranda et al, Nat. Phys. 7, 502 (2011). 

[9] D. Wang et al, Phys. Rev. Lett. 93, 243005 (2004); J. M. 
Sage et al, ibid. 94 203001 (2005); T. Rieger et al, ibid. 
95, 173002 (2005); J. Deiglmayr et al, ibid. 101, 133004 
(2008); S. D. Kraft et al, J. Phys. B 39, S993 (2006); S. 
Y. T. van de Meerakker et al, Science 4, 595 (2008); J. 
Deiglmayr et al, Phys. Rev. A 82, 032503 (2010); A. D. 
Lercher et at., Eur. J. Phys. D 1434, 1 (2011). 

[10] D.-W. Wang, M. D. Lukin and E. Dernier, Phys. Rev. 
Lett. 97, 180413 (2006). 

[II] C. Kollath, J.S. Meyer, and T. Giamarchi, Phys. Rev. 



5 



Lett. 100, 130403 (2008); A. Arguelles and L. Santos, 
Phys. Rev. A 75, 053613 (2007); C.-M. Chang et ai, ibid 

79, 053630 (2009). 

[12] A. Pikovski et al, Phys. Rev. Lett. 105, 215302 (2010); 
A. C. Potter et ai, ibid. 105, 220406 (2010); M. Baranov 
et ai, Phys. Rev. A 83, 043602 (2011); R. M. Lutchyn, E. 
Rossi and S. Das Sarma, ibid. 82, 061604 (2010)(R). 

[13] E. Burovski, G. Orso and T. Jolicoeur, Phys. Rev. Lett. 
103, 215301 (2009); G. Orso, E. Burovski and T. Joli- 
coeur, ibid. 104, 065301(2010); G. Roux, E. Burovski and 
T. Jolicoeur, Phys. Rev. A 83, 053618 (2011). 

[14] M. Dalmonte, G. Pupillo and P. Zoller, Phys. Rev. Lett. 
105, 140401 (2010). 

[15] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 

80, 885 (2008). 

[16] Coupled ID wires can be obtained, e.g., by selective evap- 
oration using (electric) field gradients, analogous to [171 ] . 
Amodsen Chotia, private communications. 

[17] S. Stock et al, Phys. Rev. Lett. 95, 190403 (2005). 

[18] S. A. Meek et ai, Phys. Rev. A 83, 033413 (2011); S. A. 
Meek et ai, Science 324, 1699 (2009). 

[19] T. Giamarchi, Quantum Physics in one dimension, (Ox- 
ford University press, Oxford, 2003). 

[20] See EPAPS Document No. [] for further details. 

[21] E. Haller et ai, Nature 466, 597 (2010). 

[22] T. Lahaye et ai, Rep. Prog. Phys. 72, 126401 (2009). 

[23] U. Schollwock, Rev. Mod. Phys. 77, 259 (2005). 

[24] D. S. Petrov et al, Phys. Rev. Lett. 99, 130407 (2007). 

[25] A. Micheli et al, Phys. Rev. A 76, 043604 (2007). 

[26] B. Friedrich and D. Herschbach, Phys. Rev. Lett. 74, 
4623 (1995); S. Kotochigova and E. Tiesinga, Phys. Rev. 
A 73, 041405(R) (2006). 

[27] C. Weitenberg et al, Nature 471, 319-324 (2011); J. Si- 
mon et ai, ibid. 472, 307-312 (2011). 

[28] B. Wunsch et al, Phys. Rev. Lett. 107, 073201 (2011). 



SUPPLEMENTARY MATERIAL 
Luttinger staircase of dimer crystals 

In the balanced, strongly interacting regime, where 
CDW correlations dominates over DL ones, the system 
may undergo a BKT phase transition towards a crys- 
talline, fully gapped phase in presence of an optical lat- 
tice. When the periodic potential lattice spacing A/2 is 
a integer multiple of the mean interparticle distance 1/n, 
that is n\/2 = 1//, / € N, the optical lattice is properly 
described by a sine-Gordon term[l|: 

H b 6h a U a n\d x d 2 J dxcos[2f ^2 Mx)l (2) 

where the coefficient U a is proportional to the depth of 
the lattice. According to BKT scaling, a crystalline phase 
is stabilized if the TLL parameter of the charge sector is 
smaller than a cr itical value K { c c) = 4// 2 [l|, thus form- 
ing a Luttinger staircase in analogy with the single tube 
case [H • Each DC is characterized by a periodic structure 
with one dimer every / sites and is fully gapped. 

Balanced case: numerical results 

In this section, we provide additional informations 
about the numerical calculation in the balanced case. 
The first transition of interest is the CDW-DC one, which 
corresponds to the pinning of composite bosons by the 
underlying lattice. In both regions, A s > 0, whereas the 
charge sector is gapless in the CDW and gapped in the 
DC. A clear signature of such a transition is given by an 
exponential decay of T>(x); however, due to the large ex- 
ponent of such a correlation close to the transition line, 
we choose to apply an alternative method based on the 
investigation of the so called block entropy (BE) : 

S(l) = -Trp A \og 2 p A (3) 

where p A is the reduced density matrix of a block A 
of length I. For critical model subject to open bound- 
ary conditions, and described by a conformal field theory 
(CFT) , it was shown that (see [3j] and references therein) : 

5^(0 = ^^(0 + const + O(l/0 (4) 
6 

where c is the so called central charge, L is the system size 
and / the cord distance. In our case, DL and CDW phases 
are described by a CFT with c = 1 (see Fig. H}, whereas 
the DC is a non-critical, fully gapped phase whose block 
entropy saturates to a constant value over a certain dis- 
tance. As proposed in [3| , the BE can be used to charac- 
terize a transition from critical to non-critical phase by 
computing 



AS(L) = S L (L/2) - S L/2 (L/4). 



(5) 



6 



In the thermodynamic limit, this quantity saturates to 
c/6 for critical theories and to for non-critical ones: red 
points and blue squares in Fig. 2 reflect this property, 
as can be inferred from the numerical data presented in 
Fig. \E\. We notice that oscillations typical of Tomonaga- 
Luttinger liquids Q do not allow for an exact estimate of 
c, even though the difference between critical and non- 
critic 




FIG. 4: (color online). Entanglement entropy of a subblock of 
length I with respect to the system for d = 0.5, 1.4: numerical 
datas (blue points and red squares respectively) and fits to 
CFT prediction which give a central charge c = 0.981, 0.97, 
in good agreement with the predicted value for a gapless phase 
phase with A s ^ 0. 




FIG. 6: (color online). Correlation functions T>(x),Q(x) 
(red circles and blue triangles respectively) in the CDW (left 
panel) and DL phase (right panel) obtained from DMRG in a 
L — 100 system. Here, g = 0.5 and d — 0.5, 1.4 respectively; 
the dashed line in the left panel denotes a reference curve 




FIG. 7: (color online). Density distribution (ni,i +712,;) in the 
balanced case for ni = 712 = 1/5; orange diamonds, blue cir- 
cles, red triangles and green squares denote d — 0.5, 1.3, 2 and 
2.3 respectively. The density distribution becomes sharper 
and sharper as we approach the DC phase. 
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FIG. 5: (color online). Scaling of AS as a function of the sys- 
tem size; points denote numerical datas, whereas thick lines 
are fits. The dashed line correspond to c = 1. For d < 2.2, the 
asymptotic is in good agreement with a c — 1 gapless phase, 
whereas it goes to zero in the DC phase. 



The DL-CDW transition is instead more subtle, as 
it corresponds to a precise value of the TLL parameter 
K c = 1. While an accurate determination of the phase 
boundary is involved, we just characterized point in the 
phase diagram by looking at T>(x) and Q{x). Sample 
numerical datas representing points close to the phase 
boundary are presented in Fig. [5] . 



Finally, it is worth mentioning that the density distri- 
bution changes drastically as we increase the interaction 
strength. In Fig. [Jj we notice that in the DC parti- 
cles are locked to certain positions, whereas in both gap- 
less phases the density shows modulations which become 
smaller in the DL phase. 



Unbalanced case 

Unbalanced gases have been shown to give rise to a se- 
ries of exotic pairing instabilities when the density ratios 
of the two component is a rational number In gen- 

eral, the TLL parameter of both components should be 
very small in order to establish such exotic bound states; 
in our particular case, one can estimate the TLL param- 
eter of the n + 1 composite field Dm = — '?2)/v / 2 
by neglecting effective spin-charge interactions, getting 
K M = { k2k \ +k - ){1 + kK 1 K 2X /{^K 1 +K 2 ))- 1 ' 2 , where 
we have considered equal velocities V\ = v 2 and a gen- 
eral forward scattering term proportional to g' = —vx- 
This results recovers previous expressions presented in 
Refs.jBI, and can be extended to consider general den- 
sity ratios n\jn 2 G Q. However, the stability of bound 
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FIG. 9: (color online). Phase diagram in the unbalanced case 
FIG. 8: (color online). Density distributions m(blue circles) n 2 = 2ni = 1/10 at fixed d 2 = 0.5: black squares and red 
and n 2 (red triangles) in the phase separated state (see text). circles denote 2TLL and TL phase respectively. The dashed 

line is a guide for the eye. 



states with more that three particles is challenging; in 
fact, large density imbalances may favor phase separa- 
tion or collapse over pairing. This intriguing perspective 
goes beyond the scope of the current work, in which only 
TLs were investigated via numerical simulations. 



Numerical Results 

We discuss now some additional numerical results on 
the specific unbalanced case n 2 = 2n\. In Fig. [5] we plot 
the density distribution of both species in the strongly 
interacting regime, d\ = 5.6,^2 = 0.9, t% = 0.1, where 
we have fixed t\ = 1 and n 2 = 1/10; DMRG calcula- 
tions were performed in this regime by taking up to 512 
states per block and 7 sweeps in L = 120 chain. The sys- 
tem phase separates into small high density islands where 
ri\ ~ n 2 and regions where ri\ = 0. This phenomenon is 
justified by the fact that such strong interaction induces 
a strong repulsion between heavy particles, whereas the 
interwire interaction is strongly attractive, thus favoring 
two-particle pairing, as present in the islands. 

Let us turn our attention on the role of the interaction 
unbalance d 2 jd\. From the analytical treatment, it is 
difficult to extract quantitative informations about how 
large this quantity should be in order to realize a trimer 
phase. In FigOJ we present the phase diagram of the sys- 
tem at a fixed interaction strength d 2 = 0.5 and different 
interaction asymmetry 1 < dijdy < 8. The TL phase 
survives even for relatively small interaction unbalance, 
albeit the critical hopping ratio £2/^1 becomes smaller 
and smaller as we approach d\ = 2d 2 . 



Density locking in the TL phase 

The choice of a local effective order parameter related 
to the TL is of great interest when the goal is to investi- 
gate an inhomogencous situation such as a trapped gas. 



In case of multiparticle pairing, a natural choice, in anal- 
ogy with the two-component caseQ, is S K n = kti\ — n 2l 
which in the trimer case becomes: 

Sn = 2ni — n 2 . (6) 

In Fig. Uni we plot Sn for different values of t 2 /ti with 
di = 3 = 6d 2 and n 2 = 2n\ = 1/10 in a L = 80 sites 
chain; the oscillatory profile, due to open boundary con- 
ditions, display a maximum smaller than Sn c ~ 0.05 
when we are in the TL phase (blue and orange lines), 
whereas it's considerably larger for t 2 > 0.2. 

A direct comparison between homogeneous and inho- 
mogeneous setup can be performed by focusing on the 
central part of the system. In Fig. [TT] (left panel), we plot 
Sn for a L = 120 system around the trap minimum: for 
large values of the trapping potential, the density lock- 
ing is not as tight as needed in order to stabilize a TL in 
the homogeneous case, whereas for more shallow traps, 
JC < 0.1, the oscillations are compatible with a TL. As a 
further check, we considered correlation function decay in 
the central region of the system, plotted in FigfTT] (right 
panel) in the shallow trap regime; C and T> correlations 
decay exponentially, only T being algebraically decaying, 
confirming the presence of a TL in the central part of the 
system. 

Coexistence of dimers and trimers 

In presence of slightly different density ratios with re- 
spect to exact commensurabilities the sine-Gordon term 
responsible for exotic pairing instabilities 0, Q does os- 
cillate, and thus bound states in general are not sta- 
ble. However, the spatial-dependent potential drastically 
changes this picture. In fact, particles can rearrange 
in order to realize a commensurate phase in the mid- 
dle of the system, expelling particles in excess towards 
the edges. 
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FIG. 10: (color online). Density locking Sn in the homoge- 
neous case, calculated in a L = 80 site chain. Green squares, 
red triangles, blue circles and orange stars represent different 
values of ia/ii = 0.3,0.2,0.1,0.05, the other relevant param- 
eters being n 2 — 2n± = 1/20, di = 3 = 6d2- 




file of both species for K, = 0.03, t% = 0.1, d\ = 3 = 6c?2 in 
a L = 160 sites system with unbalanced, non commen- 
surate densities, Ni = 9,N2 — 14. The trap center still 
presents a strong density locking of 2??i — ri2, whereas in 
the outer part the system resemble a DL where ni = n-?: 
in th 
exist 




FIG. 12: (color online). Density distribution ni (black line), 
7i2(red line) and 2*rn (blue points) for K, = 0.03, ti = 0.1, d\ = 
3 = 6^2 in a L = 160 sites system, with N\ = 9, N2 = 14; 
in the central part, n2 — 2ni, whereas away from the center 
ni ~ n2. 



FIG. 11: (color online). Panel a: density locking in the in- 
homogeneous case: K, = l(green diamonds), 0.3 (black cir- 
cles), 0.1 (red squares) and 0.03(blue triangles). Panel 6: 
correlation functions in the inhomogeneous case, JC = 0.03, 
in a L = 120 site system, with respect to x = L/4; T(x) 
(dot-dashed), T>(x) (black) and Ci(x) (red). In both panels, 
di = 3 = 6d 2 ,t 2 = 0.1. 



The 2ni < n-i case is plotted in [8[ : in this case, trimers 
form in the middle, whereas extra heavy particles occupy 
both wings. In the opposite limit 2n2 > 2n\ > ri2 one 
can expect naively a similar picture, specially because 
of the much smaller mass of light particles with respect 
to trimers. However, when intermediate-to-large interac- 
tions are considered (e?i * d 2 ^ ti), this configuration is 
energetically unfavored because of relevant interwire at- 
traction which prevents light particles to occupy region of 
space where ni — 0. In Fig. [12] we plot the density pro- 
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